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Abstract. We use Virtual Observatory (VO) tools to identify optically faint, obscured (i.e., type 2) active galactic 
nuclei (AGN) in the two Great Observatories Origins Deep Survey (GOODS) fields. By employing publicly 
available X-ray and optical data and catalogues we discover 68 type 2 AGN candidates. The X-ray powers of 
these sources are estimated by using a previously known correlation between X-ray luminosity and X-ray-to- 
optical flux ratio. Thirty-one of our candidates have high estimated powers (L x > 10 44 erg/s) and therefore 
qualify as optically obscured quasars, the so-called "QSO 2". Based on the derived X-ray powers, our candidates 
are likely to be at relatively high redshifts, z ~ 3, with the QSO 2 at z ~ 4. By going ~ 3 magnitudes fainter 
than previously known type 2 AGN in the two GOODS fields we are sampling a region of redshift - power space 
which was previously unreachable with classical methods. Our method brings to 40 the number of QSO 2 in the 
GOODS fields, an improvement of a factor ~ 4 when compared to the only 9 such sources previously known. We 
derive a QSO 2 surface density down to 10 -15 erg cm -2 s _1 in the 0.5 — 8 keV band of <; 330 deg -2 , ~ 30% of 
which is made up of previously known sources. This is larger than current estimates and some predictions and 
suggests that the surface density of QSO 2 at faint flux limits has been underestimated. This work demonstrates 
that VO tools are mature enough to produce cutting-edge science results by exploiting astronomical data beyond 
"classical" identification limits (R ^ 25) with interoperable tools for statistical identification of sources using 
multiwavelength information. 
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1. Introduction that is optically obscured, radio-quiet type 2 QSO, ex- 

ist. Indeed, until very recently, very few, if any, examples 
The unified model for active galactic nuclei (AGN) is of thig dags were knQwn Apart from thdr importance for 

largely accepted (e.g., Urry & Padovani EMI see also the AGN modelS; type 2 gources are expected to make a signif _ 

very recent results by Jaffe et al. WDM • The apparent dis- icant fraction of the x _ ray background (see; e . g . ( Comastri 

parate properties and nomenclature of active galaxies can et ^ gjjjjjj and are therefore also cosmologically very rele- 

be explained by the physics of black hole, accretion disk, yant _ Thege gources are heavily reddened and there fore fall 

jet, and obscuring torus convolved with the geometry of through the « stand ard" (optical) methods of quasar selec- 

the viewing angle. Type 1 sources are those in which we tion ^ Thc hard x _ rayg; however; are thought to be able to 

have an unimpeded view of the central regions and there- penetrate the t0 rus. Type 2 QS0; therefore, should have 

fore exhibit the straight physics of AGN with no absorp- nam)W; if any ^ permit t e d lines (and might look like nor- 

tion. Type 2 objects arise when the view is obscured by mal galaxies in the optical/UV band), powerful hard X-ray 

the torus. While many examples of local, and therefore emissiorii andj in some caseSj a high equivalent width Fc 

relatively low-power, type 2 AGN are known (the Seyfert R Hne ^ Norman et aL JM)^ . 
2s), it has been debated if their high-power counterparts, 

In this paper we use Virtual Observatory (VO) 

Send offprint requests to: Paolo Padovani tools to identify 68 type 2 AGN candidates in the two 
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Great Observatories Origins Deep Survey (GOODS) fields 
(Giavalisco et al. 2004a), ~ 1/2 of which qualify as 
QSO 2 candidates. Based on the properties of already 
known sources, we expect the large majority of these to 
be obscured quasars whose identification is only possible 
through their X-ray emission. 

VO initiatives are now at a stage where prototype tools 
can be utilised to produce scientific results. Real gains 
have been made in the areas of accessing and describing 
remote data sets, manipulating image and catalogue data, 
and performing remote calculations in a fashion similar 
to grid computing. These prototype tools are enabled by 
the VO infrastructure and interoperability standards that 
are being developed cooperatively by all the VO projects 
under the auspices of the IVOA 1 . VO software is expected 
to mature significantly over the next 1-2 years as the VO 
projects progress from demonstrations to building robust 
systems. In this paper we have taken advantage of the first 
interoperability gains to produce the "first science" for the 
VO. 

In Section 2 we present the data we have used, Section 
3 discusses the tools we employed, while Section 4 de- 
scribes our method. Section 5 presents our results, which 
are discussed in Section 6, while Section 7 summarizes 
our conclusions. Throughout this paper spectral indices 
are written S v oc v~ a and we adopt a cosmological model 
with H = 70 km s" 1 Mpc" 1 , Q M = 0.3, and fi A = 0.7. 

2. The data 

We use for our purposes the two GOODS fields (Giavalisco 
et al. l2004af> . namely the Hubble Deep Field-North (HDF- 
N) and the Chandra Deep Field-South (CDF-S). Since 
GOODS includes some of the deepest observations from 
space- and ground-based facilities, these are the most 
data-rich, deep survey areas on the sky. The GOODS 
field centres (J2000.0) are 12 ft 36 m 55 8 , +62°14'15" for the 
HDF-N and 3 h 32 m 30 s , -27°48'20" for the CDF-S. Each 
field provides an area of approximately 10' x 16'. 

Deep X-ray (Chandra) catalogues are available for a 
larger region around both fields (Alexander et al. 120031 
Giacconi et al. 2002). For consistency, we use the cata- 
logues produced by Alexander et al. |2003) for both the 
2 Ms HDF-N and 1 Ms CDF-S data. These include 503 
(HDF-N) and 326 (CDF-S) objects respectively, for a to- 
tal of 829 sources. The data cover the 0.5 — 8.0 keV band 
and the catalogues provide counts and fluxes in various 
sub-bands. Note that, due to the twice as long exposure 
time, the HDF-N reaches fainter X-ray fluxes and includes 
a (54%) larger number of sources. 

In the optical we use the publicly available GOODS 
Hubble Space Telescope (HST) Advanced Camera for 
Surveys (ACS) data (proposal ID 9 583). The observations 
consist of imaging in the F435W, F606W, F775W, and 
F850LP passbands, hereafter referred to as B, V, i, and z, 
respectively. We use here version vl.O of the reduced, cal- 



ibrated, stacked, and mosaiced images and catalogues as 
made available by the GOODS team 2 . 

Finally, identifications and redshifts through optical 
spectroscopy are available from Barger et al. IJ2003II and 
Szokoly et al. <PU4) for the HDF-N and CDF-S respec- 
tively. We note that ~ 56% of the X-ray sources in 
the GOODS fields have spectroscopic redshift determina- 
tions. These sources are necessarily relatively bright, with 
(I) ~ 21 (HDF-N; only 3% fainter than 24th mag) and 
(R) ~ 22 (CDF-S; only 5% fainter than 25th mag). 

3. Virtual Observatory Tools 

Astronomy is at a turning point. Major breakthroughs 
in telescope, detector, and computer technology allow as- 
tronomical surveys to produce massive amounts of im- 
ages, spectra, and catalogues. These datasets cover the sky 
at all wavelengths from 7- and X-rays, optical, infrared, 
through to radio. The VO is an international, community- 
based initiative, to allow global electronic access to avail- 
able astronomical data, both space- and ground-based. 
The VO aims also to enable data analysis techniques 
through a coordinating entity that will provide common 
standards, wide-network bandwidth, and state-of-the-art 
analysis tools. 

The Astrophysical Virtual Observatory Project 
(AVO) 3 is conducting a research and demonstration pro- 
gramme on the scientific requirements and technolo- 
gies necessary to build a VO for European astronomy. 
The AVO has been jointly funded by the European 
Commission (under FP5 - Fifth Framework Programme) 
with six European organisations participating in a three 
year Phase-A work programme. 

The AVO project is driven by its strategy of regular 
scientific demonstrations of VO technology. For this pur- 
pose an "AVO prototype" has been built. The prototype 
consists of a suite of interoperable software, plus a set of 
conventions or standards for accessing remote data, and 
for launching remote calculations. The main component 
of the software is based on the CDS Aladin visualisation 
interface (Bonnarel et al. 2000). This prototype VO por- 
tal (v. 1.003-/3) allows efficient interactive manipulation of 
image and catalogue data, and provides access to remote 
data archives and image servers via the GLU registry of 
services 4 . 

The Aladin image server is an example of such 
a VO service. It describes the images stored in the 
Aladin database using a data model (Images Distributees 
Heterogenes pour l'Astronomie; IDHA 5 ), and provides im- 
age cutouts on request. In this paper we make heavy use 
of cutouts of the GOODS data available via this service. 
The prototype is also interoperable with the other long 
standing CDS Vizier and SIMBAD services (Ochsenbein, 
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Bauer, & Marcout l2l)l)l)[) . and significant interoperability 
gains are achieved by use of the VOTable 6 format for as- 
tronomical tables. 

The prototype includes a catalogue cross matching ser- 
vice. This service allows positional cross matching of two 
catalogues to find the best matched source, all matching 
sources, or sources not matching within a given threshold 
radius (in arcseconds). These three modes, plus the ability 
to compare directly with the images from which the cata- 
logues were generated, make for an extremely efficient tool 
for which to perform cross matches and check for multiple, 
or aberrant matches. 

In addition to the AVO software, we have also made 
intensive use of the Star link topcat tool , and the VOIndia 
VOPlot plugin 8 to the prototype. 



4. The method 

The two key physical properties that we use to identify 
type 2 AGN candidates are that they be obscured, and 
that they have sufficiently high power to be classed as an 
AGN and not a starburst. To find these candidates, we use 
a relatively simple method based on the X-ray and optical 
fluxes. 

Optical data alone are not sufficient for this purpose 
because at the redshifts that we expect to find these 
sources the nuclear, AGN emission can be diluted by 
the host galaxy. Indeed, Moran, Filippenko & Chornock 
( 2002 ) have shown that ~ 60% of local Seyfert 2 galaxies 
would be classified as normal galaxies if no decomposition 
of their optical emission were available. However, AGN re- 
veal themselves by their hard X-ray emission and power. 
Ideally, to find type 2 AGN, one would need X-ray spec- 
tra to select sources with flat spectral indices, indicative of 
absorption. The typical count rates of the X-ray sources in 
the GOODS fields are however relatively low for detailed 
spectral fitting and therefore we select sources based on 
their X-ray hardness ratio, a measure of the fraction of 
hard photons relative to soft photons. 

To estimate the X-ray power for candidate type 2 
sources we use the correlation described by Fiore et al. 
PTOljl between the f (2 - WkeV) / f (R) ratio and the hard 
X-ray power (see their Fig. 5). The basis of this correla- 
tion is the fact that the f(2 — 10keV)/f(R) ratio is roughly 
equivalent to the ratio between the nuclear X-ray power 
and the host galaxy R band luminosity. Since the host 
galaxy R band luminosities (unlike the X-ray power) show 
only a modest amount of scatter, this flux ratio is a good 
indicator of X-ray power. 

Each of the steps in our method, including our tech- 
nique for identifying the optical counterparts of the X-ray 
sources, are described in detail below. 
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Fig. 1. The distribution of hardness ratios in the 
Alexander et al. (2003) catalogues. We define as absorbed 
sources to the right of the dashed line, that is having hard- 
ness ratio HR > —0.2. 



4.1. Selecting absorbed sources 
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The Alexander et al. I|2003J) catalogues provide counts 
in various X-ray bands. We define the hardness ratio 
HR = (H - S)/{H + S), where H is the hard X-ray 
counts (2.0 — 8.0 keV) and S is the soft X-ray counts 
(0.5 - 2.0 keV). Following Szokoly et al. PM|l when H 
is an upper limit we set HR = — 1, while when S is an 
upper limit we set HR = 1. The 14 sources with upper 
limits in both H and S bands were excluded. Szokoly et 
al. I|2004[) have shown that absorbed, type 2 AGN are 
characterized by HR > —0.2. We adopt this criterion 
and identify those sources which have HR > —0.2 as ab- 
sorbed sources. We find 294 (CDF-S: 104, HDF-N: 190) 
such absorbed sources which represent 35^2% (lc Poisson 
errors are from Gehrels 1986) of the X-ray sources in the 
Alexander catalogues. The hardness ratio distribution is 
shown in Fig. ^ We note that for an a x — 0.8 the flux in 
the 2 — 8 keV band used by Alexander et al. ( 2003 ) is only 
~ 15% smaller than that in the commonly used 2 — 10 keV 
band. We also point out that increasing redshift makes the 
sources softer (e.g., at z — 3 the rest-frame 2 — 8 keV band 
shifts to 0.5 — 2 keV) so our selection criterion will mis- 
takenly discard some high-z type 2 sources, as pointed out 
by Szokoly et al. (2004). The number of type 2 candidates 
we find has therefore to be considered a lower limit. 

One commonly adopted definition of absorbed source 
is TVh > 10 22 cm~ 2 , where iVn is the intrinsic absorption 
at the redshift of the source. We discuss below how this 
compares with our definition. 
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4.2. Finding the optical counterparts 

The optical counterparts to the X-ray sources were se- 
lected by cross-matching the absorbed X-ray sources with 
the GOODS ACS catalogues (29 599 sources in the CDF- 
S, 32 048 in the HDF-N). The GOODS catalogues contain 
sources that were detected in the z-band, with BVi pho- 
tometry in matched apertures (Giavalisco et al. [2004b). 
We checked for possible offsets between the optical and 
X-ray astrometry by cross-matching the full (CDF-S and 
HDF-N) Alexander catalogues with the GOODS cata- 
logues using a threshold radius of 1.25" (see below). In 
the HDF-N we detected a systematic shift of —0.029" in 
R.A., and —0.297" in declination of the Alexander posi- 
tions with respect to the GOODS positions, and we ap- 
plied this correction to the Alexander coordinates before 
doing the cross-match. 

Note that the ACS image areas for the CDF-S and 
the HDF-N are both smaller than (and are completely 
within) their respective Chandra image fields. 546 out 
of the 829 sources in the Alexander catalogue (CDF-S: 
222/326, HDF-N: 324/503) fall within the ACS image 
area. Of the 294 absorbed sources, 203 (CDF-S: 77, HDF- 
N: 126) fall within the ACS image area. 

To find the optical counterparts of the absorbed X- 
ray sources, we initially search for optical sources that lie 
within a relatively large threshold radius of 3.5" (corre- 
sponding to the maximal 3a positional uncertainty of the 
X-ray positions) around each X-ray source. This was done 
using the cross match facility in the AVO prototype tool 
using the "best match" mode. 

We find 195 "best" matches, 76 in the CDF-S and 119 
in the HDF-N. So, in the CDF-S, all but one of the ab- 
sorbed sources in the ACS image area have optical coun- 
terparts. The remaining source, Alexander ID 213, falls 
within the 10" disk of a bright galaxy. In the HDF-N there 
are there are 7 absorbed sources within the image area 
that do not have an optical counterpart (these 7 sources 
do however have optical matches in Barger et al. l|2"003). 
although we note that most of these are below the 5<r limit 
of their catalogue). 

Since the 3.5" radius is large relative to the median po- 
sitional error, and given the optical source density the ini- 
tial cross match inevitably includes a number of false and 
multiple matches. To limit our sample to good matches, we 
use the criterion that the cross match distance be less than 
the combined optical and X-ray 3<r positional uncertainty 
for each individual match. Applying this distance /error 
< 1 criterion we limit the number of matches to 168 (CDF- 
S: 65, HDF-N: 103). These matches are all within a much 
smaller radius than our initial 3.5" threshold, with most 
of the distance/error < 1 matches being within 1.25" (and 
two matches at 1.4 and 1.5"). 

Considering not only the best match but also the mul- 
tiple matches within a threshold radius of 1.25" (and dis- 
tance/error < 1) we find 189 matches (CDF-S: 67, HDF- 
N: 122) to the 168 X-ray sources. This means that our 



method of only considering the best matches does discard 
21 possibly valid matches in preference for a closer match. 

The number of false matches we expect to have, as 
detailed in Appendix lAl is small, between 8 and 15%. 

Our method requires R magnitudes in order to esti- 
mate the X-ray luminosity via the Fiore et al. correlation. 
The ACS band closest to the R band is the i band. We 
then convert the ACS i magnitudes to the R band assum- 
ing (R — iacs) = 0.5, which is the typical value we derive 
for previously known type 2 AGN in the GOODS fields 
(see Sect. I6.1f) . The R band flux, f(R), was then com- 
puted by converting R magnitudes to specific fluxes and 
then multiplying by the width of the R filter (Zombeck 
1990), as in Fiore et al. (HEH. 

4.3. Estimating the X-ray power 

As discussed above, previously classified sources and their 
spectroscopic redshifts are available from two catalogues: 
Szokoly et al. $M)ty for the CDF-S and Barger et al. 
(2003) for the HDF-N. For these sources we derived the 
2 — 8 keV X-ray power, L 2 -$, using our adopted cosmology 
and assuming cv x — 0.8 for the k-correction (Fiore et al. 
2003). No absorption correction was applied, which means 
that the intrinsic powers of our sources could be larger. As 
we are dealing with the hard X-ray band, however, this is 
probably not going to make much difference, as pointed 
out by Fiore et al. ffflty . 

The known sources in the HDF-N were easily identi- 
fied in our list of candidates because Barger et al. ( 2003 ) 
used the same Alexander et al. (2003) catalogue of X- 
ray sources. All 103 HDF-N sources that passed the dis- 
tance/error < 1 criterion are listed in Barger et al. (2003), 
54 of which have spectroscopic redshifts, leaving 47 unclas- 
sified HDF-N sources. For the CDF-S Szokoly et al. piHljl 
utilised the X-ray catalogue of Giacconi et al. (2002) 
[which is derived from the same data as Alexander et al. 
( 2003 ) but with somewhat different source extraction pro- 
cedures] , and list only the positions of their optical coun- 
terparts as determined from VLT FORS1 R-band imaging. 
To establish which CDF-S candidates are already listed in 
the Szokoly et al. (]2()()4p catalogue we cross matched the 
GOODS ACS optical positions of our candidates with the 
Szokoly optical positions, and inspected all cases where 
there were multiple or possibly confused matches. Out of 
the 65 CDF-S sources which passed the distance/error < 1 
criterion, 54 are listed in Szokoly et al. (20(2}> and 44 of 
these have spectroscopic redshifts. This leaves a total of 
68 (CDF-S: 21, HDF-N: 47) unclassified candidates. 

For the unclassified sources we estimated the X-ray 
power as follows: we first derived the /(2 — lOkeV) / f (R) 
flux ratio, and then estimated the X-ray power from 
the correlation found by Fiore et al. (2003), namely 
logL 2 -io = logf{2 - lOkeV)/ f(R) + 43.05 (Fiore, p.c; 
see their Fig. 5). Note that this correlation has an r.m.s. 
of ~ 0.5 dex in X-ray power and that, since the X-ray 
powers in the Fiore et al. I|2003|) correlation have been 
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corrected for absorption, the estimated powers are already 
automatically corrected. 

We stress that our estimated X-ray powers reach ~ 
10 45 erg/s and therefore fall within the range of the Fiore 
et al. C 2003J) correlation. 

4.4. Finding the type 2 AGN candidates 



The work of Szokoly et al. ( 2004) has shown that absorbed, 
type 2 AGN are characterized by HR > —0.2. It is also 
well known that normal galaxies, irrespective of their mor- 
phology, have X-ray powers that reach, at most, L x < 
10 42 erg/s (e.g., Forman, Jones & Tucker IT9941 Shapley, 
Fabbiano & Eskridge IMTTl Cohen I2UU3J 1 . Therefore, any 
X-ray source with HR > —0.2 and L x > 10 42 erg/s should 
be an obscured AGN. Furthermore, following Szokoly et 
al. i|2004jl . any such source having L x > 10 44 erg/s will 
qualify as a type 2 QSO. 

It is a well-established fact that, at a given optical 
magnitude, galaxies are also weaker X-ray emitters than 
AGN, that is they have much smaller X-ray-to-optical flux 
ratios, typically ^0.1 (Maccacaro et al. H988fl . Indeed, the 
Fiore et al. I|2003[) correlation implies that L x < 10 42 erg/s 
corresponds to /(2 — 10keV)/f(R) < 0.1, so a selection 
based on X-ray-to-optical flux ratio would produce the 
same results. 



5. Results 

Our selection criteria have produced 68 type 2 AGN candi- 
dates, based on the definitions and the method described 
above. These sources are listed in Tab. ^ (HDF-N) and 
H (CDF-S). Tabledgives the GOODS ID in column (1), 
the Alexander et al. (2003) IDs in column (2), the X-ray 
position in columns (3)- (4), the offset between the X-ray 
and optical positions in arcseconds in column (5), the full- 
band (0.5 — 8 keV) X-ray flux in column (6), the hardness 
ratio in column (7), the ACS i magnitude in column (8), 
the estimated X-ray power in column (9). Tab. |21 gives, 
in addition, the Szokoly al. I|2()()4jl IDs in column (2), and 
the remaining columns are then shifted by one. 

Out of our candidates, 31 satisfy the further require- 
ment Z/2-io > 10 44 erg/s and therefore qualify as QSO 
2 candidates. We note that the distribution of estimated 
X-ray power covers the range 5 x 10 42 — 2 x 10 45 erg/s 
and peaks around 10 44 erg/s (see Fig. Q3 dashed line). 
The number of QSO 2 candidates, therefore, is very sensi- 
tive to the dividing line between low- and high-luminosity 
AGN, which is clearly arbitrary and cosmology depen- 
dent. For example, if one defines as QSO 2 all sources with 
Li-io > 5 x 10 43 erg/s, a value only a factor of 2 below the 
commonly used one and corresponding to the break in the 
AGN X-ray luminosity function (Norman et al. l2002|) . the 
number of such sources increases by ~ 50%. We also note 
that, based on the r.m.s. around the Fiore et al. (2003) 
correlation, the number of QSO 2 candidates fluctuates in 
the 13 — 54 region. The number of type 2 AGN, on the 



other hand, can only increase, as all our candidates have 
estimated log £2-10 > 42.5. 

It is interesting that all candidates with HR > —0.2 
have estimated X-ray power L x>eat > 10 42 erg/s, and 
therefore all qualify as AGN. Some previously known 
sources, however, do have X-ray powers below this value 
(see Sect. lB~T|) . 

As expected, being still unidentified, our sources are 
very faint: their median ACS i magnitude is ~ 25.5, which 
corresponds to R ~ 26 (compare this to the R ~ 22 typ- 
ical of the CDF-S sources with redshift determination). 
The QSO 2 candidates are even fainter, with median i 
magnitude ~ 26.3 (R ~ 26.8). 

It is important to notice that 6 of the 21 sources in Tab. 
El have actually been observed by Szokoly et al. (2004; sec 
their Fig. 5) but either their spectrum is featureless (4 
sources) or the continuum is so weak that no extraction 
could be made (2 sources). 

5.1. Testing our estimated X-ray powers 

To check the reliability of our estimated X-ray powers we 
compared the predicted and observed luminosities for the 
already known sources. The Fiore et al. (20Q2J) correla- 
tion is valid for non-type 1 AGN and galaxies so we did 
the comparison for these sources. Namely, we excluded 
the type 1 sources from the Szokoly et al. l|2*0"04) cat- 
alogue and the sources with broad lines (type "B") in 
Barger et al. (2003|. Moreover, we only considered sources 
with quality flag > 2.0 in the former catalogue and ex- 
cluded objects with less secure redshifts (type "s") in 
the latter, to have reliable redshifts and therefore powers. 
This leaves us with 122 objects. For these, we find that 
£x,est <x £ x 93±0 - 06 , where the slope is the mean ordinary 
least-squares slope (Isobe et al. H990|) . The mean values of 
the estimated ((log£ Xj6S t) = 42.57 ± 0.08) and observed 
((logL x ) = 42.49 ±0.09) luminosities, are consistent. This 
shows that our method works for non-type 1 AGN and 
can then be safely applied to optically unidentified but 
absorbed sources. 

It is tempting to use our estimated X-ray powers to- 
gether with the observed fluxes to derive redshifts for our 
type 2 candidates. We can first check how reliable these 
are by comparing them with the spectroscopic redshifts 
available for the known non-type 1 AGN described above. 
We follow standard procedures (see, e.g., Mobasher et al. 
2004) and quantify the reliability of the estimated red- 
shifts, z es t, by measuring the fractional error for each 
galaxy, A = (z ost - z spC c)(l + z spoc ). We find a median 
error, (A) = 0.005, an r.m.s. scatter, <r(A) ~ 0.29, and a 
rate of "catastrophic" outliers, 77, defined as the fraction 
of the full sample that has |A| > 0.2, ~ 19%. Clipping 
these outliers gives (A) = —0.03 and er(A) ~ 0.1. For 
comparison, Mobasher et al. Ij2004(l have used extensive 
multiwavelength photometric data to estimate photomet- 
ric redshifts for a sample of 434 galaxies with spectro- 
scopic redshifts in the CDF-S. They find (A) = -0.01, 
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Table 1. Type 2 AGN Candidates, HDF-N. 



GOODS ID 



A03 



RA 



DEC 



OFFSET 
arcsec 



/(0.5 

io- 1_ 



SkeV) 
erg cm -2 s~ 



HR 



log £ x (2 - SkeV) 
erg/s 



J123542.91+621144.6 
J123556. 10+621219. 6 
J123606.32+621233.2 
J123606.86+621021.7 
J123611. 78+621015.0 
J123613. 73+620800.5 
J123614.46+621045.5 
J123616.05+621108.1 
J123616.12+621514.0 
J123620.88+621415.8 
J123621.96+621603.6 
J123622. 59+621306.5 
J123622.64+621630.2 
J123624.81+621744.2 
J123627.60+621158.8 
J123634.46+620942.3 
J123634.62+621422.1 
J123635.28+621152.2 
J123637.12+620636.7 
J123644.13+621245.2 
J123646.58+620857.0 
J123647.94+621020.6 
J123651.28+621052.0 
J123653. 17+621947.8 
J123654.11+620835.4 
J123654.59+621111.8 
J123658. 77+621022.7 
J123700.45+621509.2 
J123702.44+621926.0 
J123706.88+622225.6 
J123707.21+621408.5 
J123709.93+622259.4 
J123711.34+621330.9 
J123712.00+621325.9 
J123712.04+621212.3 
J123713.84+621458.7 
J123714.68+621839.7 
J123720.34+621353.6 
J123723.99+621304.7 
J123726.51+622026.8 
J123726.85+621747.9 
J123731. 13+621030.1 
J123731. 54+621305.9 
J123734.09+621624.8 
J123737.06+621834.5 
J123737.95+621713.9 
J123804.23+621526.0 



24 


12 35 42.98 


+62 11 44.4 


0.50 


48 


12 35 56.14 


+62 12 19.2 


0.33 


75 


12 36 06.38 


+62 12 32.2 


0.84 


79 


12 36 06.84 


+62 10 21.4 


0.11 


92 


12 36 11.80 


+62 10 14.5 


0.28 


99 


12 36 13.73 


+62 08 00.2 


0.03 


102 


12 36 14.45 


+62 10 45.0 


0.15 


108 


12 36 16.03 


+62 11 07.7 


0.18 


109 


12 36 16.11 


+62 15 13.7 


0.01 


123 


12 36 20.91 


+62 14 15.5 


0.20 


129 


12 36 21.94 


+62 16 03.8 


0.48 


133 


12 36 22.57 


+62 13 06.2 


0.10 


135 


12 36 22.66 


+62 16 29.8 


0.23 


143 


12 36 24.78 


+62 17 43.7 


0.32 


152 


12 36 27.75 


+62 11 58.4 


1.10 


181 


12 36 34.48 


+62 09 41.8 


0.32 


185 


12 36 34.64 


+62 14 22.1 


0.32 


188 


12 36 35.26 


+62 11 51.6 


0.32 


198 


12 36 37.26 


+62 06 37.5 


1.48 


228 


12 36 44.12 


+62 12 44.5 


0.39 


242 


12 36 46.50 


+62 08 56.9 


0.59 


246 


12 36 47.94 


+62 10 19.9 


0.41 


266 


12 36 51.28 


+62 10 51.5 


0.20 


276 


12 36 53.07 


+62 19 47.8 


0.69 


281 


12 36 54.13 


+62 08 35.1 


0.17 


284 


12 36 54.58 


+62 11 10.6 


0.88 


303 


12 36 58.79 


+62 10 22.2 


0.30 


312 


12 37 00.45 


+62 15 08.8 


0.08 


321 


12 37 02.43 


+62 19 26.1 


0.44 


341 


12 37 06.72 


+62 22 24.5 


1.35 


347 


12 37 07.20 


+62 14 07.9 


0.30 


362 


12 37 09.88 


+62 22 58.7 


0.50 


366 


12 37 11.38 


+62 13 30.8 


0.38 


368 


12 37 12.04 


+62 13 25.7 


0.35 


369 


12 37 12.09 


+62 12 11.3 


0.84 


375 


12 37 13.90 


+62 14 58.0 


0.60 


380 


12 37 14.60 


+62 18 39.5 


0.57 


402 


12 37 20.34 


+62 13 53.1 


0.20 


412 


12 37 24.00 


+62 13 04.3 


0.13 


423 


12 37 26.51 


+62 20 26.8 


0.27 


424 


12 37 26.77 


+62 17 47.8 


0.56 


430 


12 37 31.13 


+62 10 30.0 


0.22 


431 


12 37 31.55 


+62 13 06.1 


0.54 


434 


12 37 33.98 


+62 16 24.0 


0.93 


445 


12 37 37.04 


+62 18 34.4 


0.24 


447 


12 37 37.90 


+62 17 13.3 


0.42 


486 


12 38 04.11 


+62 15 25.2 


1.00 



1.51 


0.23 


27.51 ± 0.20 


44.8 


18.7 


0.78 


22.50 ± 0.02 


43.9 


0.96 


0.14 


25.58 ± 0.07 


43.8 


0.74 


1.00 


24.56 ± 0.07 


43.4 


2.50 


0.69 


25.94 ± 0.18 


44.5 


1.62 


0.09 


25.54 ± 0.17 


44.0 


3.63 


0.54 


24.70 ± 0.06 


44.1 


16.07 


0.27 


25.16 ± 0.10 


44.9 


1.02 


-0.09 


26.22 ± 0.15 


44.1 


0.96 


0.46 


24.76 ± 0.08 


43.5 


1.22 


1.00 


27.33 ± 0.27 


44.7 


0.28 


-0.13 


24.70 ± 0.04 


42.9 


1.02 


1.00 


24.66 ± 0.07 


43.6 


0.47 


1.00 


25.78 ± 0.14 


43.7 


2.22 


0.50 


23.87 ± 0.03 


43.6 


0.57 


-0.10 


27.12 ± 0.20 


44.1 


0.12 


0.33 


26.23 ± 0.14 


43.3 


0.40 


-0.07 


24.66 ± 0.05 


43.0 


0.25 


1.00 


27.01 ± 0.17 


43.9 


0.08 


1.00 


24.99 ± 0.10 


42.9 


0.51 


1.00 


23.75 ± 0.04 


43.0 


1.93 


0.54 


26.83 ± 0.27 


44.7 


2.42 


0.77 


24.06 ± 0.04 


43.7 


0.64 


1.00 


27.00 ± 0.18 


44.4 


0.66 


1.00 


22.50 ± 0.01 


42.6 


0.94 


0.20 


26.55 ± 0.10 


44.2 


1.56 


0.69 


24.88 ± 0.07 


43.8 


1.80 


-0.19 


25.17 ± 0.08 


43.9 


3.97 


-0.18 


26.93 ± 0.24 


44.9 


0.39 


1.00 


25.22 ± 0.09 


43.3 


0.98 


0.23 


25.39 ± 0.11 


43.8 


5.21 


0.64 


23.90 ± 0.07 


43.9 


0.56 


0.11 


24.38 ± 0.07 


43.1 


0.92 


1.00 


25.11 ± 0.10 


43.7 


0.38 


0.33 


25.94 ± 0.14 


43.6 


1.11 


1.00 


26.50 ± 0.19 


44.4 


0.48 


1.00 


25.62 ± 0.06 


43.6 


1.35 


1.00 


23.02 ± 0.02 


43.0 


6.07 


0.49 


24.62 ± 0.04 


44.3 


1.40 


0.17 


24.99 ± 0.07 


43.8 


1.68 


0.17 


26.10 ± 0.18 


44.3 


0.54 


-0.06 


25.31 ± 0.10 


43.4 


0.40 


-0.18 


26.04 ± 0.13 


43.6 


0.82 


1.00 


25.42 ± 0.08 


43.8 


4.43 


-0.13 


27.18 ± 0.29 


45.1 


0.86 


0.22 


24.78 ± 0.10 


43.5 


2.77 


0.43 


22.58 ± 0.03 


43.2 
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Table 2. Type 2 AGN Candidates, HDF-S. 



GOODS ID 


S04 


A03 


RA 


DEC 


OFFSET 
arcsec 


/(0.5- 
10~ 15 erg 


■ 8keV) 

— 2 — I 

cm s 


HR 


i 


logL x (2-8fey) 
erg/s 


J033205.35-274644.1 




50 


03 32 05.35 


-27 46 44.1 


0.04 




3.37 


0.58 


26.11 ± 0.19 


44.7 


J033206. 12-274927.8 


259 


52 


03 32 06.13 


-27 49 27.7 


0.16 




4.43 


0.64 


23.81 ± 0.05 


43.8 


J033214.86-275103.8 


82a 


100 


03 32 14.87 


-27 51 03.8 


0.17 




1.74 


0.02 


25.82 ± 0.16 


44.2 


J033217.11-274136.9 




116 


03 32 17.11 


-27 41 36.6 


0.28 




1.60 


0.18 


26.32 ± 0.15 


44.4 


J033219.95-274243.1 


524 


130 


03 32 19.96 


-27 42 43.2 


0.10 




1.43 


0.03 


27.18 ± 0.20 


44.7 


J033221. 78-274656.5 




138 


03 32 21.75 


-27 46 57.0 


0.66 




0.15 


1.00 


26.63 ± 0.24 


43.7 


J033222. 54-274603.8 


145 


145 


03 32 22.54 


-27 46 03.9 


0.06 




4.74 


0.32 


24.59 ± 0.04 


44.2 


J033224.54-274010.4 




156 


03 32 24.55 


-27 40 10.7 


0.28 




0.89 


1.00 


22.54 ± 0.01 


42.8 


J033224.85-275600.0 




159 


03 32 24.84 


-27 56 00.0 


0.08 




5.61 


-0.05 


23.74 ± 0.04 


43.9 


J033226. 78-274604.0 


518 


175 


03 32 26.75 


-27 46 04.3 


0.49 




0.74 


0.11 


27.44 ± 0.35 


44.5 


J033229. 78-275328.4 




187 


03 32 29.84 


-27 53 28.6 


0.78 




0.77 


1.00 


27.03 ± 0.23 


44.4 


J033232. 16-274651.5 


515 


201 


03 32 32.17 


-27 46 51.4 


0.10 




1.57 


0.39 


27.01 ± 0.24 


44.7 


J033233. 28-274236.0 




206 


03 32 33.27 


-27 42 36.1 


0.21 




4.01 


0.43 


25.51 ± 0.11 


44.5 


J033235. 21-275318.0 


148 


218 


03 32 35.23 


-27 53 17.8 


0.37 




3.19 


0.08 


27.10 ± 0.27 


45.0 


J033236. 16-275037.1 


C09 


223 


03 32 36.17 


-27 50 37.2 


0.22 




0.68 


1.00 


25.40 ± 0.11 


43.8 


J033238. 76-275121.6 




237 


03 32 38.78 


-27 51 21.6 


0.25 




0.92 


0.12 


26.15 ± 0.08 


44.1 


J033239.05-274439.3 




241 


03 32 39.06 


-27 44 39.1 


0.25 




2.64 


-0.11 


26.64 ± 0.19 


44.7 


J033239. 17-274832.4 




243 


03 32 39.17 


-27 48 32.7 


0.33 




0.20 


1.00 


26.27 ± 0.13 


43.8 


J033243.04-274845.1 


256 


255 


03 32 43.05 


-27 48 45.0 


0.19 




4.05 


0.69 


25.15 ± 0.06 


44.3 


J033246.34-274632.0 


147 


268 


03 32 46.35 


-27 46 31.9 


0.13 




7.03 


0.64 


24.95 ± 0.11 


44.5 


J033247.06-275333.5 




272 


03 32 47.05 


-27 53 33.3 


0.26 




3.02 


-0.03 


26.21 ± 0.13 


44.6 



<r(A) ~ 0.1, and r\ ~ 2.4 — 10%, depending on the sub- 
sample considered. Excluding the outliers they find values 
of cr(A) ~ 0.05. Given the simplicity of our method, these 
results are certainly encouraging and show that our esti- 
mates should at least provide a rough idea of the redshifts 
at which we might expect our sources to be. 

Applying this method to our type 2 AGN candidate list 
we find a mean of (z cs t) ~ 2.9 (median 2.6). The QSO 2 
have, as expected, significantly higher estimated redshifts 
(zest) ~ 3.7 compared to (z cst ) ~ 2.2 for the non-QSO 
candidates. 

We can also use the ACS four band photometry to 
constrain the redshifts of our sources. Cristiani et al. 



(2004) have discussed a set of four optical criteria us- 
ing the i — z, B — V, and V — i colours (a variation of 
the "B-dropout" technique) to select AGN in the redshift 
range 3.5 Ss z Ss 5.2. By applying the same criteria to 
our candidates we find 9 such sources. Their average es- 
timated redshift is (z cst ) ~ 3.6, and all but one sources 
have 2.8 < z cst < 5.4. Again, this shows that our esti- 
mated redshifts are, overall, relatively robust. 

These sources are listed in Table which gives the 
GOODS ID in column (1), the ACS i magnitude in column 
(2), followed by the i — z, V — i and B — V colours (all in 
the AB system) in columns (3)- (5). 
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HERE FIG. 2 

Fig. 2. CDF-S Cutouts. B, V, i, and z ACS image cutouts are displayed for all the type 2 AGN candidate in the CDF-S. 
Each image is 3"x 3"in size, with north up and east to the left. The Alexander identification number is shown in the 
top left of each panel, and the circles drawn on the i-band images indicate the Alexander X-ray positions, with the 
radius indicating the 90% positional error. The (black and white) square symbols indicate the positions of the GOODS 
source which have been matched to the X-ray positions. The X symbols indicate the positions of any corresponding 
sources in the Szokoly catalogue. 

HERE FIG. 3 

Fig. 3. CDFN Cutouts. B, V, i, and z ACS image cutouts are displayed for all the type 2 AGN candidate in the CDFN. 
Each image is 3"x 3"in size, with north up and east to the left. The Alexander identification number is shown in the 
top left of each panel, and the circles drawn on the z-band images indicates the Alexander X-ray positions, with the 
radius indicating the 90% positional error. The (black and white) square symbols indicate the positions of the GOODS 
sources which have been matched to the X-ray positions. 



Table 3. B dropouts. 



GOODS ID 


i 


i — z 


V -i 


B-V 


J033226.78-274604.0 


27.44 


0.37 


3.21 


-1.36 


J033232. 16-274651.5 


27.01 


0.27 


0.43 


1.25 


J033238.76-275121.6 


26.15 


0.04 


0.23 


1.41 


J033239.05-274439.3 


26.64 


0.26 


0.23 


4.41 


J123611. 78+621015.0 


25.94 


0.71 


2.46 


>3.0 


J123627.59+621158.8 


23.86 


0.24 


0.69 


1.25 


J123634.46+620942.3 


27.12 


0.28 


1.53 


>3.0 


J123714.68+621839.7 


25.62 


0.13 


0.44 


1.87 


J123731. 54+621305.9 


26.04 


-0.13 


0.46 


1.39 



5.2. Testing our method 

Out of the 9 QSO 2 previously known (see Sect. 0), 
our method has rediscovered them all as type 2 AGN, 
8/9 with L x . cs t > 8 x 10 43 erg/s (the ninth one having 
£ x ,est ~ 4 x 10 42 erg/s). Out of the 29 type 2 AGN with 
L2-8 > 10 42 erg/s in the Szokoly et al. (2004) catalogue, 
28 have L x ,cst > 10 42 erg/s and therefore fulfil our se- 
lection criteria, while the last one is barely below with 
£x,est - 7 x 10 41 erg/s. 

During the completion of this work three redshift 
databases covering the CDF-S and HDF-N became pub- 
lic. Namely: the ESO/GOODS FORS2 spectroscopy data 
(Vanzella et al. I2DU1P ). the VIMOS VLT Deep Survey 
(VVDS; Le Fevre et al. HEHH and the Team Keck 
Treasury Redshift Survey (TKRS; Wirth et al. I2TMR . 
We cross-matched our candidate list with these surveys 
using the AVO prototype, finding only one match, namely 
GOODSJ123556. 10+621219.6, with a quite featureless 
spectrum and redshift z = 0.9585 (our estimated redshift 
based on its X-ray power and flux is z — 0.93). Given the 
limits of the three surveys, that is zab < 24.5, Jab < 24, 
and R < 24.4 respectively, this is not surprising. The only 
match we found, in fact, is the brightest of our northern 
candidates. 



http : //archive . eso . org/wdb/wdb/vo/goods/f orm 

http : //cencosw . oamp . f r/ 

http : //www2 . keck . hawaii . edu/science/tksurvey/ 



5.3. Image cutouts 

Image cutouts of the ACS BViz data centred on each of 
the new AGN type 2 candidates were generated using the 
AVO prototype. This was done by sending requests to the 
Aladin image server which remotely extracts the required 
image sections in the four bands from a database con- 
taining the GOODS data. The hierarchical description of 
the ACS images in terms of the IDHA data model, and 
advanced protocols for sending detailed requests to the 
server makes this process more efficient and flexible than 
traditional image servers. Figs. [21 and |21 display 3"x3" sec- 
tions of the cutouts in all four bands for the CDF-S and 
HDF-N respectively. Each image is shown with a linear 
grey-scale over the range of minimum to maximum im- 
age values of the individual 3"x3" image sections. The 
90% error circle of the Alexander et al. 1)2003 ) catalogue 
X-ray positions is overlaid on the i band images (since 
our method used the i band), along with the GOODS 
catalogue position (black/ white squares). The Szokoly et 
al. I)2()()4|l optical positions (X symbols) are also shown 
for the CDF-S sources (Fig. EJ. Note that in a few cases 
the position of the GOODS optical counterpart is outside 
the X-ray 90% error circle. These are still valid matches 
in terms of the distance/error parameter (see Sect. 14. 4|) . 
which limits matches to within the 3ct total positional un- 
certainty which is ~ 1.8 times larger than the 90% radius. 
Our sources cover a large spectrum of morphologies, 
ranging from extended, low-surface brightness to point- 
like. It is interesting to note, however, that all QSO 2 
candidates are point-like. 



5.4. Hubble Ultra Deep Field data 

A number of sources also fall within the Hubble Ultra 
Deep Field (UDF; Beckwith et al. l2TM}l . The UDF pro- 
vides extremely deep ACS Wide Field Channel BViz 
imaging (proposal IDs 9 978 and 10 086) of a single 5.25' x 
5.25' field centred on 03 /l 32 m 39.0 s , -27°47'29.1" (J2000) 
which lies within the GOODS CDF-S field. Five of the 
absorbed sources fall within this field, and 3 of these (A03 
201, 243 and 255) are in our list of unclassified type 2 AGN 
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Table 4. Type 2 AGN candidates, UDF 



S04 A03 UDF ID 



RA 



DEC 



OFFSET 



logL x (2-8fceF) 
erg/s 



515 



256 



201 7326 
243 1441 
255 1025 



03 32 32.17 
03 32 39.19 
03 32 43.04 



-27 46 51.49 0.10 
-27 48 32.87 0.31 
-27 48 45.08 0.18 



26.869 ± 0.021 
28.328 ± 0.046 
25.153 ± 0.006 



44.6 
44.7 
44.3 



201 B 


V 


© 


z 










243 B 


V 


© 


z 










255 B 

* 
■ 


V 

« 

• 
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• 



Fig. 4. CDF-S UDF Cutouts. B, V, i, and z ACS image 
cutouts are displayed for the type 2 AGN candidates that 
fall within the area of the UDF. The (black and white) 
square symbols indicate the positions of the UDF sources 
in the UDF catalogue. The X symbols indicate the posi- 
tions of any corresponding sources in the Szokoly cata- 
logue. 

candidates (Tab. [5J. For these sources we cross matched 
the X-ray positions with the UDF Version 1.0 BViz i- 
band detected catalogues 12 to find the UDF optical coun- 
terparts. We then used the UDF i band magnitudes to 
estimate the X-ray power using the Fiore et al. (2003) 
correlation as described above. The estimated X-ray pow- 
ers, and UDF i band magnitudes for these 3 sources are 
listed in Tabled and UDF image cutouts overlaid with the 
UDF, Alexander et al. lf2TM|) . and Szokoly et al. pM)l 
positions are shown in Fig. 0J Note that for one source, 
A03 243, we find a different, fainter, optical counterpart 
than was found using the GOODS catalogue, resulting 
in a larger estimated X-ray power. Classically, one would 
have solved this dilemma by obtaining a spectrum of both 
sources. However, given the faintness of the two candi- 
dates (i ~ 26.3 and 28.3), this is unlikely to happen in the 
near future. [The two previously classified sources that fall 
within the UDF are A03 245 (QSO-2) and 242 (AGN-1)]. 
Note that none of our candidates have spectra obtained 
with the ACS High Resolution Channel grism 13 as part of 
the UDF observing campaign. 

6. Discussion 

6.1. Previously known type 2 AGN 

We defined a sample of previously known type 2 AGN 
in the GOODS ACS fields as follows. In the HDF-N, we 



1 http://www.stsci.edu/hst/udf 

1 http : //www . stecf . org/UDF/goodsdata . html 



selected absorbed (HR > —0.2) sources without broad 
lines (type "B" in Barger et al. l2003p having X-ray power 
L 2 - S > 10 42 erg/s (see Sect. EM- In the CDF-S we ini- 
tially selected sources with "AGN-2" or "QSO-2" clas- 
sifications in Szokoly et al. I|2004[) . We then applied the 
further criteria that L2-8 > 10 42 erg/s and HR > —0.2 in 
order to be consistent with the HDF-N sample. Successive 
power and HR cuts exclude 2 and 3 sources respec- 
tively, with a further 3 sources excluded because they have 
no counterpart in the Alexander et al. (2003) catalogue 
(which is necessary for consistent definitions of HR and 
-^2-s)- [Note that differences between our HR values and 
those listed by Szokoly et al. H2004fl can be explained by 
the different input catalogues produced by Giacconi et al. 
PUU2)l and Alexander et al. Q2003[l ]. These definitions lead 
to a total of 79 known type 2 AGN (CDF-S: 35, HDF-N: 
44) . These are listed in Tables and El which give the 
Alexander et al. (fSn77^1 and Szokoly et al. pTMl (for the 
CDF-S) IDs, redshift, hard X-ray power, and notes. Only 
9 of these sources can be classified as QSO 2 (£2-8 > 10 44 
erg/s) , 3 of which have poor redshift determinations (two 
sources with type "s" flags in Barger et al. 120031 and one 
source with quality flag = 1 in Szokoly et al. l2004|l . We also 
note that with a lower X-ray power cut of £2-8 > 5 x 10 43 
erg/s the number of QSO 2 increases to 13. 

Fig. shows the X-ray power distribution for our new 
type 2 AGN candidates (dashed line), previously known 
type 2 AGN (solid line), and the combined sample (dot- 
ted line) . It is interesting to note how the distributions are 
very different, with the already known type 2 AGN peak- 
ing around L x ~ 10 43 erg/s and declining for luminosities 
above ~3x 10 43 erg/s, while our new candidates are ris- 
ing in this range and peak around L x ~ 10 44 erg/s. To 
be more quantitative, while only ~ 1/5 of already known 
type 2 AGN have logL x > 43.5, ~ 3/4 of our candidates 
are above this value. This difference is easily explained 
by our use of the X-ray-to-optical flux ratios to estimate 
X-ray powers and by the fact that our candidates are on 
average ~ 3 magnitudes fainter than previously known 
sources. Our method is then filling a gap in the luminosity 
distribution, which becomes almost constant in the range 
10 42 ^ Lj ^ 3 x 10 44 erg/s. This also explains the fact 
that the number of QSO 2 candidates we find is ^ 3 times 
larger than the previously known ones. 

Conversely the fluxes of the already known QSO 2 
reach /(2 — 8keV) ~3x 10~ 15 erg cm~ 2 s _1 , while our 
new QSO 2 candidates go down to the fainter limit of 
/(2 - 8keV) ~4x 10~ 16 erg cm~ 2 s _1 . This is explained 
by the much fainter optical magnitudes we are probing 
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Table 5. Previously Known Type 2 AGN, HDF-N. Table 6. Previously Known Type 2 AGN, CDF-S. 
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42.27 




122 


1.338 


42.90 


s 


127 


1.014 


43.27 


m 


142 


0.747 


42.47 




150 


0.762 


42.82 




157 


1.264 


43.89 




158 


1.013 


43.02 




160 


0.847 


42.90 


m 


163 


0.485 


42.72 




164 


0.953 
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42.27 




190 


2.005 
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42.23 


s 


240 


0.961 


43.95 




259 


1.609 


44.48 


m 


261 


0.902 
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43.32 
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Barger 


et al. |27I03l 


s, less secure redshift identii 



A03 S04 redshift 



based primarily on a single emission line and the continuum 
shape; m, complex or multiple structure or possible contami- 
nation of the photometry by a nearby brighter object. 



with our method which, for a given X-ray-to-optical flux 
ratio, translate into fainter X-ray fluxes as well. 

By going fainter we are also probing the type 2 AGN 
population at higher redshifts. While our estimated red- 
shift is (z cst ) ~ 2.9 (median 2.6), previously known type 
2 AGN have (z) ~ 1.1 (median 0.9). For QSO 2 we find 



logL x (2-8fceV) 
erg/s 



Notes' 1 
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"From Szokoly et al. (2004): QF — 1, spectrum cannot be iden- 
tified securely, typically only a single narrow line present; QF 
= 0.5, only a hint of some spectral feature. 



(z C st) ~ 3.7 (median 3.5), as compared to (z) ~ 2.3 (me- 
dian 2.6). 

6.2. Total number of type 2 AGN 

We find a total of 147 type 2 AGN in the GOODS ACS 
area, 79 of which were already known. This corresponds 
to 27^2% of the 546 sources in the Alexander et al. (2003) 
catalogue (14^ 2 % including the previously known sources 
only). As regards QSO 2, the total number is 40, only 9 of 
which were previously known. Our method has therefore 
more than quadrupled the number of such sources in the 
GOODS ACS fields. These represent 7±{% of all (546) X- 
ray sources (the previously known QSO 2 make up 2~^\%). 
Note that all of these sources satisfy the commonly 
adopted definition of absorbed source discussed above. In 
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Fig. 5. The X-ray power distribution for our new type 
2 AGN candidates (dashed line), previously known type 
2 AGN (solid line), and the sum of the two populations 
(dotted line). QSO 2 are defined, somewhat arbitrarily, as 
having L 2 -iokcV > 10 44 erg/s. 



fact, Fig. 12 of Szokoly et al. (2004) shows that HR > 
—0.2 corresponds to Ah > 10 22 cm -2 , for an intrinsic 
a x = 1, for z <; 0.4. The lowest redshift for previously 
known type 2 AGN is 0.4, while the lowest estimated red- 
shift for our type 2 candidates is 0.9. The typical HR and 
redshifts for our sources imply that we are dealing with 
Ah ~ 10 23 and w 3 x 10 23 cm -2 for the known and can- 
didate sources respectively. 

X-ray obscuration can be also present in a small frac- 
tion of broad-lined AGN. Perola et al. (20QU an( A refer- 
ences therein) find this to be the case in ~ 10% of their 
sources, with estimated Ah <; 10 22 cm" 2 . As discussed in 
the previous paragraph, these sources are then likely to 
be at relatively low redshift. Therefore, contamination by 
type 1 AGN in our sample is expected to be negligible. 

6.3. The surface density of QSO 2 

The detection of faint type 2 AGN candidates and the 
careful assessment of the previously known such sources 
in the two GOODS fields allow us to put strong constraints 
on the surface density of type 2 sources, and in particu- 
lar on that of QSO 2. We use a value for the Viz-ba,nd 
coverage of 365 arcmin 2 ~ 0.1014 deg 2 (Giavalisco et al. 

rcnraa . 

Due to the variable sensitivity of the Chandra ACIS-I 
detector across the field of view, the area in which fainter 
X-ray sources could be detected is smaller than that for 
brighter sources. In other words, the sky coverage is not 



uniform at all fluxes. Fig. 5 of Giacconi et al. (|20(J2B and 
Fig. 19 of Alexander et al. (2003J), however, show that the 
effect is strong only at relatively small fluxes. For example, 
for both CDF-S and HDF-N the effective area decreases in 
the hard band, for which this effect is strongest, by <; 25% 
only for /(2 - 8/10keV) £ 10 -15 erg cm" 2 s _1 . 

To better quantify the magnitude of this correction 
we evaluated the sky coverage for the combined GOODS 
fields by simply scaling the two sky coverages so that the 
maximum area for each field was equal to half the to- 
tal GOODS area, summing then up the two areas. This 
works only as a first approximation and slightly overesti- 
mates the correction since the GOODS ACS fields are at 
the centre of the Chandra fields and therefore deeper than 
average in the X-ray band. In other words, the sky cov- 
erage should be re-computed as its shape would change, 
with a larger area accessible at fainter fluxes and therefore 
a smaller correction than we estimate. In any case, un- 
der our assumptions we find that the correction is within 
the lcr Poisson range, and therefore within the statistical 
uncertainties, for /(0.5 — 8keV) ^ 10 -15 erg cm -2 s -1 . 
Our conservative approach is then to not to take the sky 
coverage into account and limit ourselves to X-ray fluxes 
> 10 -15 erg cm -2 s -1 . This means that our values are 
actually lower limits, although the real surface densities 
should not be more than 25% larger. 

Table [7| gives the QSO 2 surface density for four flux 
limits, splitting the contribution into known and candidate 
sources for the full 0.5 — 8 keV band. (The counts in the 
hard band are not very different, since we are dealing with 
absorbed sources, with the total numbers changing to 33 
[from 34] and 8 [from 10] for the first and third flux limit 
respectively.) 

These numbers can be compared with recent estimates 
and predictions. Perola et al. (J2004) find a surface den- 
sity of highly obscured QSO, which they define as having 
L 2 _io > 10 44 erg/s and A H > 10 22 cm" 2 , ~ 48 deg -2 for 
/(2 — lOkeV) > 10 -14 erg cm -2 s -1 , consistent with our 
value. At lower fluxes the situation is different. For exam- 
ple, for f(0.5-7keV) > 5 x 10 -15 erg cm -2 s -1 Gandhi et 
al. (2004) quote an estimated value > 3 deg -2 and possi- 
bly higher than 10 deg -2 and a predicted value from their 
model of 19 deg -2 . We find a density ~ 100 deg" 2 , that 
is a factor > 10 larger than their estimate and ~ 5 times 
larger than their prediction. Already known sources make 
up ~ 70% of our total value, so this high density is clearly 
not attributable only to our new candidates. We do note 
that 5 out of the 7 previously known QSO 2 come from 
the paper by Szokoly et al. (2004|, which appeared only 
recently and could not be accounted for by Gandhi et al 
(2004). For the same flux limit and L2-10 > 3 x !° 44 er g/ s 
these authors predict a surface density ~ 3 deg -2 , while 
we find 5 sources (3 of them already known), which im- 
plies a density 49^2? deg -2 . The decrease in space density 
at high luminosity is therefore much less than predicted 
by their model, namely only a factor ~ 2 instead of ~ 6. 

At even fainter fluxes, Gandhi et al. (2004) quote a 
value from their model of 35 deg -2 above /(0.5 — 7keV) ~ 
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Table 7. QSO 2 surface density. 



/(0.5 - 8keV) 



erg cm 2 s 1 



1 x 10 



2 x KT 15 
5 x 10~ 15 
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138^36 227^47 

30+f 6 99+^ 
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25 

14 

3 

1 



34 

23 



89 
89 
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4 30 



2 x 10 15 erg cm 2 s 1 , while we find a density ~ 230 
deg -2 , that is a factor ~ 6 larger than their prediction. 

The population synthesis model of Gilli, Salvati & 
Hasinger <|2001l and private communication) predicts 16 
obscured QSO with intrinsic luminosity L2-10 > 3 x 10 44 
erg/s, N H > 10 22 cm" 2 , and z > 3 in the 1 Ms CDF-S. 
We find 2 known plus 7 candidate sources above these red- 
shift and power limits with HR > —0.2. According to Fig. 
12 of Szokoly et al. IJ2TMJ1 at z > 3 an A H = 10 22 cm" 2 
corresponds to HR ~ —0.6 for an intrinsic a x — 1, which 
means that our definition is more restrictive. Considering 
also that the southern GOODS ACS area is ~ 60% smaller 
than the CDF-S area, there are strong indications that our 
number might be a very solid lower limit to such sources 
and that therefore our findings are not inconsistent with 
the numbers predicted by Gilli et al. H2001fl . 

As mentioned in Sect.EJ the number of our QSO 2 can- 
didates depends on the log L2-10 ~ l°g /(2 — 10keV) /f(R) 
correlation of Fiore et al. (2003), which has an r.m.s. of 
~ 0.5 dex in X-ray power. However, even if we consider 
the worst case scenario in which all our estimated X-ray 
powers are too large by 0.5 dex the number of QSO 2 
candidates for the four flux limits in Tab. are 13, 8, 2, 
and 1 respectively. In other words, our total densities de- 
crease, at worst, by ~ 25 — 35% at the lowest fluxes, and 
are basically unchanged at larger fluxes. 

The resolved fraction of the X-ray background due to 
QSO 2 down to /(2 - 8keV) = 10~ 15 erg cm" 2 s" 1 is 
estimated to be 10+ 2 %. We have used here the value of the 
total X-ray background measured by UHURU and HEAO- 
1 (Marshall et al. I198U|I and the uncertainties reflect the 
r.m.s. of the Fiore et al. (2003) correlation (the Poissonian 
error is even smaller) . Note that given the relatively small 
area of the GOODS fields we are missing the bright end of 
the number counts. In fact, only four of our sources have 
/(2 — 8keV) > 10 -14 erg cm -2 s" 1 and these reach only 
2.6 x 10 -14 erg cm~ 2 s _1 . This fact, the discussion above 
and the points detailed in the next section, all mean that 
our value has to be regarded as a robust lower limit. 



6.4. Caveats and comments 

In this paper we have employed a statistical method to 
identify very faint type 2 AGN. We had also to rely on an 
empirical technique to estimate the X-ray powers, which 
were needed to make sure that the sources we are dealing 



with have AGN-like outputs. As such, we can only pro- 
vide a list of candidates and not firm classifications. This 
was obviously expected, as the great majority of our candi- 
dates are so faint that even the largest telescopes presently 
available would require extremely long exposures to secure 
a decent spectrum. 

However, we believe that our method is robust, 
as shown by the various checks we have performed. 
Importantly, we have been very conservative in our es- 
timates of the number of type 2 AGN candidates and the 
surface densities we have estimated need to be considered 
lower limits (modulo what discussed in the previous para- 
graph), for the following four reasons: 1. our selection cri- 
terion for absorbed sources (HR > —0.2) will mistakenly 
discard some high-z type 2 sources f Sect. EH} : 2. the same 
criterion is also more restrictive than the commonly used 
one based on Nh fSect. l6~2"|) : 3. the X-ray powers for pre- 
viously known sources have not been corrected for absorp- 
tion, which means more sources could be above the QSO 
2 limit (Sect- IO)! ; and finally, 4. sky coverage effects, once 
properly taken into account, will reduce the available area 
at smaller X-ray fluxes and therefore increase the source 
surface density fSect. l6~3|) . 

6.5. Next steps enabled with VO tools 

As noted above, a remaining source of uncertainty with 
the newly discovered sample is the reliance on an empirical 
technique in the determination of the X-ray power of the 
objects. 

In Sect. 15.11 we noted how the 4 colour ACS photo- 
metric data have been used to estimate redshifts for a 
number of objects in a restricted range. Future work will 
enable the use of the 4 band ACS plus IR (VLT/ISAAC 
and Spitzer) photometry to determine photometric red- 
shifts over the full GOODS field. VO technologies are 
to be employed to facilitate this work. Upgrades to the 
AVO capability will include the provision of a "redshift- 
determination" service, which will automate the gener- 
ation of point spread function matched multi-band in- 
put photometric catalogues via the use of AstroGrid's 
(Walton, Lawrence & Linde I2UU3JI ACE (Astronomical 
Catalogue Extractor) service [a Web service providing 
access to the SExtractor application (Bertin & Arnouts 
1996)], and feed these into a range of photometric redshift 
determination techniques [e.g., Bayesian prior method, 
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Bpz (Benitcz 2000), the neural network technique, ANNz 
(Collister & Lahav l2"oT>!}) and Hyperz (Bolzonella, Miralles 
& Pello l2lIUU|) ]. 

It is necessary to include the IR data in order to obtain 
reliable photometric redshifts, and with the additional 
photometric bands, the effects of possible dust extinction 
will also be identifiable (Rowan-Robinson I2003[l . 



7. Conclusions 

We have used Virtual Observatory tools to identify ob- 
scured AGN much fainter than previously known ones, 
going beyond what is currently possible using the "clas- 
sical" approach of classifying sources by taking spectra 
of them even using the largest telescopes available today. 
The fact that we have obtained scientifically useful and 
cutting-edge results is proof that VO tools have evolved 
beyond the demonstration level to become respectable re- 
search tools. The VO is already enabling astronomers to 
reach into new areas of parameter space with relatively 
little effort. 

Our main results can be summarized as follows: 

1. By employing publicly available, high-quality X-ray 
and optical data we have discovered 68 type 2 AGN 
candidates in the two GOODS fields. Their X-ray pow- 
ers have been estimated by using a previously known 
correlation between X-ray luminosity and X-ray-to- 
optical flux ratio. Thirty-one of our candidates have 
high luminosity (£2-iokcV > 10 44 erg/s) and therefore 
qualify as QSO 2, that is optically obscured quasars. 
Based on their derived X-ray powers, our candidates 
are likely to be at relatively high redshifts, z ~ 3, with 
the QSO 2 at z ~ 4. 

2. We have tested our method and results extensively 
and find them to be very robust. In particular: a) 
our approach recovers most (~ 97%) of the previously 
known type 2 AGN in the GOODS fields; b) the X- 
ray power estimates agree very well with those derived 
from spectroscopic redshifts for the non-type 1 AGN in 
the GOODS fields; c) the redshifts derived from our es- 
timated powers are consistent with those measured for 
previously known non-type 1 AGN; d) the spectrum of 
our brightest northern candidate, recently made pub- 
lic as part of the Team Keck Treasury Redshift Survey, 
confirms our type 2 classification as it is basically fea- 
tureless, with a measured redshift extremely close to 
our estimate. 

3. By going ~ 3 magnitudes fainter than previously 
known type 2 AGN we are sampling a region of red- 
shift - power space so far unreachable with classical 
methods. Our method brings to 40 the number of QSO 
2 in the GOODS fields, an improvement of a factor 
~ 4 when compared to the 9 such sources previously 
known. The relatively large QSO 2 number density we 
derive Q> 330 deg~ 2 for /(0.5 - 8keV) > 10~ 15 erg 
cm~ 2 s -1 , ~ 30% of which is made up of already known 



sources) suggests that its value at faint flux limits has 
been underestimated. 

For the handful of our candidates with R ^ 26 spec- 
troscopy with a class 8 - 10 m telescope is still within 
reach, albeit with relatively long (sS 10 h) exposure times. 
We plan to follow these up and confirm (or refute) their 
classification in the near future. 

In closing we note that this paper represents the first 
significant published science result that has been fully en- 
abled via end-to-end use of VO tools and systems. 

VO initiatives are in their early stages. Significant 
progress is being made both by the AVO and other na- 
tional VO projects. Construction of a truly pervasive sys- 
tem is beginning to provide the end user access to a pow- 
erful combination of data access, applications embedded 
within user workflows, running over high speed networks 
on powerful compute resources. Early use, exploiting a 
prototype system in the USA, constructed by the US-VO 
project 14 to simplify the cross matching of Two Micron 
All Sky Survey and Sloan Digital Sky Survey source cat- 
alogues has enabled Berriman et al. Il2003jl to discover a 
small number of hitherto undiscovered Brown Dwarf can- 
didates. Their work demonstrated possibilities, but did 
not alter the scientific understanding of that particular 
problem. 

In this paper we have shown how, even with the AVO's 
early prototype, enough capability and access to data is 
available to make the scientific mining of the, in this case 
GOODS, data significantly easier for the "general" as- 
tronomer. This is because gains in interoperability sim- 
plify the acquisition of the data, and enable quick inter- 
operation of a number of common, but computationally 
intensive, tasks such as cross matching and multi layer 
visualisation of both image and catalogue data. 

With the rapid evolution of the VO, science discovery 
will be routinely performed using VO techniques, an early 
example of which is described here. 

The AVO prototype used in this paper is freely avail- 
able for download 15 . 
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Appendix A: False Match Estimate 

We estimated the number of false matches using the for- 
malism described in Derriere (200Q. The expected total 

1 http://www.us-vo.org 

15 http : //www . euro-vo . org/twiki/bin/view/Avo/SwgDownload 
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distance (arcsec) 

Fig. A.l. The histogram of distances between closest 
match optical and X-ray sources overlaid with the best 
fitting model 

number of matches D, for N± objects with positional un- 
certainty a against a uniformly distributed set of objects 
with density A is: 



N 1 ^(d,X)+N^(d,X) 



D{d,N 1 ,X,a,N) 

-.Act 2 f_^_ 

2ct 2 



27rAcr 2 



■ exp 



- 1 



(A.l) 



where d is the distance between matches, $(<i, A) is 
the 2-d Poisson density function 16 and N is the number of 
objects with a true counterpart. D may also be expressed a 
sum of three components, a 17 the number of true matches, 
(3 1S the objects with a true partner but not incorrectly 
assigned and tp 19 the number of objects that do not have 
a counterpart but have been matched. 

Using our known values of N\: the number of ab- 
sorbed sources (203), A: the optical source density of the 
GOODS catalogues (0.0469 arcsec -2 ), and ct: the median 
positional error of the X-ray positions of the absorbed 
sources (0.31") we vary N to perform a maximum like- 
lihood fit of the expected distribution of matches D to 
the observed histogram of distances between closest match 
optical and X-ray sources out to the initial 3.5" threshold 
radius. The best fit occurs for N — 160. Fig. IA.ll shows 
the histogram of distances between closest match optical 
and X-ray sources overlaid with the best fitting model. 

Integrating D and a out to 1.25" we find that 92% 
of the matches out to this radius are expected to be true 
according to the model, making our false match fraction 
8%. Note that this approach neglects clustering, which 
would slightly increase this value. 

We also estimated the number of false matches by com- 
paring the number of X-ray to optical matches we find 

16 $(d, A) = 27rAdexp(-7rAd 2 ) 

17 a=5&**(d,A)exp' ■'" 



= N$(d, A) exp (- 

ip = (N!- N)$(d, A) 



using the correct coordinates, to the number of matches 
we find when a shift (of 20" much larger than the typical 
match radius and the spacing between optical sources) is 
applied to the X-ray coordinates. Within 1.25" we find 188 
matches and 38 matches when the coordinates are shifted. 
Our selection of real matches however, also included the 
distance/error criterion which means that a number of 
matching sources within 1.25" of the X-ray source were 
discarded. This criterion should also be taken into con- 
sideration when estimating the rate of false matches. We 
estimate this effect on the false match rate by calculating 
the rate of optical sources within 1.25" of the real X-ray 
position which have distance/error > 1 (10), and then 
subtract this from the rate of false matches we find in 
the offset cross-match. Therefore the false match fraction 
using this method should be (38-10)/188 = 15%. 

The fraction of false matches should then be in the 
range 8 — 15%. 
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